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Abstract. The mobility of charge carriers in a semiconductor nanowire is explored as 
a function of increasing radius, assuming low temperatures where impurity scattering 
dominates. The competition between increased cross-section and the concurrent 
increase in available scattering channels causes strongly non-monotonic dependence 
of the mobility on the radius. The mter-band scattering causes sharp declines in 
the mobility at the wire radii at which each new channel becomes available. At 
intermediate radii with the number of channels unchanged the mobility is seen to 
maintain an exponential growth even with multiple channels. We also compare the 
effects of changing the radial scaling of the impurity distribution. We use transverse 
carrier wavefunctions that are consistent with boundary conditions and demonstrate 
that the (5-function approximate transverse profile leads to errors in the case of remote 
impurities. 



1. Introduction 

The relentless miniaturization of devices has made the physics of lower dimensions a 
commonplace experimental reality today. A variety of materials are routinely fabricated 
into nanometer scale elements by research groups around the world; at low temperatures, 
nanoscale wires can physically behave as one- dimensional conductors in the sense that 
the carriers are confined to a single, or at most a limited number, of quantum modes or 
channels in the transverse directions. The study of transport in such one-dimensional 
conductors continues to reveal novel physical behavior. In this context semiconducting 
nanowires in particular have been the focus of sustained study because of the ever- 
growing technological impact of semiconductors. Several semiconducting materials are 
being employed to grow nanowires, and to fabricate elementary devices, including silicon 
(H 12] , gallium arsenide [3 E] , germanium [3 IH] , indium phosphide [71 |H] and indium 
arsenide [Oj. 

The carrier mobility is obviously a critical parameter characterizing the transport 
properties and ultimately the applications potential of a semiconductor nanowire. 
Experimental studies frequently cite the observed mobility range for samples of interest. 
For quasi one- dimensional wires made of a specific semiconducting material, the mobility 
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of the carriers depend on several factors, such as the available scattering mechanisms 
for the carriers, the carrier concentration, the temperature and the physical dimensions 
of the wire. Numerous theoretical papers [12111111121113111111311111131111 have studied 
the influence of essentially every relevant physical parameter on the mobility in ultra 
small semiconductor wires. However such studies have generally been in regimes where 
either (i) the wires are thick enough and the temperatures sufficiently high that the 
transport has bulk behavior, or (ii) they are in the extreme size quantized limit with 
only one or two channels available. Therefore it is of interest to see how the mobility in a 
semiconductor nanowire behaves when the number of available channels increases as the 
wire radius becomes larger and the system deviates from the extreme size quantization 
limit. That is our goal in this paper. A strong practical motivation exists because 
with current fabrication methods it is already possible to grow nanowires with specific 
diameters with a precision of within 10% standard deviation 5312312]) ^"^^ that is 
certain to improve rapidly with time. 

In some earlier descriptions of mobility in wires at the size quantized regime 
artificial transverse profiles were sometimes used for the carrier concentration, chosen 
primarily for mathematical convenience fTJ E| • We treat the carrier transverse motion 
quantitatively using carrier density profiles that are consistent with physical boundary 
conditions and considering multiple distinct transverse modes. Size quantization is most 
pronounced in the regime of very thin wires and low temperatures. Phonon scattering 
is suppressed at sufficiently low temperatures, so in this paper we work in a regime in 
which the mobility is mainly determined by impurity scattering. Some recent studies 
[T7j have considered impurities concentrated on the axis. Here, we instead take the 
impurities to be distributed outside the region of high carrier density since this captures 
modulation doping j22j or surface roughening effects resulting from some nanowire 
fabrication methods. 

The rest of the paper is arranged as follows: In section we develop the relaxation 
time description of the mobility applied to one dimensional systems. We specialize to 
nanowires with a circular cross-section in section 13 and derive the impurity scattering 
matrix elements. In section El we compute the radial scaling of those matrix elements 
and of the carrier mobility assuming a surface distribution of impurities. In section I3 
we present the results of our numerical calculations and our observations regarding the 
behavior of the mobility and related properties as a function of wire radius. 

2. Mobility in one dimension in terms of relaxation times 

The size quantization limit is most easily achieved at low temperatures when most 
of the carriers have insufficient energy to populate higher channels. The effect of 
phonon scattering is negligible at sufficiently low temperatures and impurity scattering 
is dominant. In order to focus exclusively on the impurity scattering mechanism we 
consider the degenerate limit corresponding to T = 0. Our description will be based on 
the relaxation time approximation to the Boltzmann equation that approximates the 
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collision term by the quotient of the deviation from equilibrium of the Fermi distribution 
function and a characteristic relaxation time. In quasi one dimension, the distribution 
function /j(k) carries a continuous label k for the wavevector along the length of the wire 
and a set of discrete channel indices, denoted by that label the degrees of freedom 
in the restricted transverse dimensions; the same labelling is used for other physical 
parameters. Solving the Boltzmann equation in the relaxation time approximation 
results in a set of equations for the carrier relaxation times Tj of each channel [221 ^ 
specific energy E 



cos(0) . (1) 



Here ^^^(k, k') is the scattering amplitude from a state in channel i with wavevector 
k to another in channel j with wavevector k'. The scattering mechanism satisfies 
J (k, k') = ^^^^(k', k). The angle (p is between the initial k and the final k' wavevectors. 

2.1. Scattering Matrix in One Dimension 

In one-dimension the wavevectors are scalars so that k = ±k allowing only two possible 
values for the scattering angle, = 0, vr. The total energy of a carrier in the ith channel 
has two parts 

E=^{k' + ^). ^^(^) = |^^'- (2) 

The longitudinal energy ei{k) is assumed parabolic with the lattice potential of the 
semiconductor incorporated into the effective mass m taken to be approximately the 
same for all channels. The transverse energy h'^K^/2m belongs to a discrete spectrum 
of energies determined by the boundary conditions of the transverse profile of the wire. 

We carry out the sum over final momenta in Eq. using the Fermi's Golden rule 
expression for the scattering amplitude 

Y.S.A^,k') = dk,\{k\j\V%i)\' 6 + - ^{k'^ + ^])) 

= T.^\{±k\J\V\k,^)\^ (3) 

Energy conservation determines the final wavevector given the initial one and the specific 
channels involved, k' = Jk"^ + [nf — k^). 



2.2. Mobility in One Dimension 

A knowledge of the scattering matrix determines the relaxation times which in turn 
determine the electron mobility for each channel in the wire through the expression 
_ (lT.kvl[ri{E)d,^(k)ft{k)] 

where the carrier velocity for the parabolic case is Vk = hk/m. The derivative is with 
respect to the energy of the free carriers. With Ep denoting the one- dimensional Fermi 
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energy, we can define Fermi energies for individual cfiannels Ei = Ep — ^f^f and tlie 
associated Fermi wavevectors. Note tliat we leave out tfie label 'F' for Fermi surface 
quantities for individual channels as superfluous since, as we see presently, we will only 
be working with Fermi surface values for each channel. At low temperatures where 
impurity scattering is most pronounced, the Fermi distribution function is essentially at 
the degenerate limit 

Mk) = e{E, - t,{k)), d,^^k)Mk) = -5{E, - e,{k)) (5) 

Allowing for spin degeneracy, the linear density of carriers in the ith channel is given 
by 

/oo die '?lc- 
^9{E, - e,{k)) = 9{E,)^ (6) 
-oo ZTT IT 

For a specific carrier density, n = N/V, and cross-section A (assumed uniform) of the 
wire the linear density is ni = nA. Then the one-dimensional Fermi energy Ep can be 
determined by adding together the densities (jHl) of each channel 

|n, = j:OiEp - ^ni)^-^Ep-K^ (7) 

i 

where 9{x) is the Heaviside unit step function. The mobility in each channel becomes 



2ki J-oo 2ti rn^ \ 2m J 

--HE.f-^ (8) 
The average mobility for electrons is then given by 

/. = 5i^ = ^(E.)^5iMM (9) 

In the degenerate limit the relaxation times have to be evaluated only at the Fermi 
surface determined by the total carrier density in the wire and the wire radius. 



3. Cylindrical nanowires 

We will now specifically consider nanowires with an uniform cylindrical cross-section, 
and a Coulomb scattering potential arising from ionized impurities. At low temperatures 
where impurity scattering is dominant a Coulomb potential is an appropriate choice 
for the scattering potential V in typical low dimensional systems. The natural basis 
functions are a product of transverse functions involving Bessel functions and a plane 
wave corresponding to the longitudinal part. So the matrix elements due to an impurity 
at ro = {zo,ro,9o} is 

{k,,j\Vir,ro)\h,t) = ^ l/drj^'de $*(r, e)$,(r, ^) 

dz (10) 
(r^ + rl- 2rro cos(6')) + {z - zqY 
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where the azimuthal angle 6 is measured from the direction of the impurity at ^q- 
Because of the cyhndrical symmetry the value of will not influence the scattering 
probability. The quantity e is the dielectric constant. It is well known ^2] that 
due to screening of the Coulomb interaction, in the strict degenerate limit of zero 
temperature the static dielectric function in one dimension evaluated in a random phase 
approximation (RPA) has a divergence at twice the Fermi vector, 2fcj, for a channel. 
From our analysis above it is clear that for intra-channel scattering at the degenerate 
limit the momentum change corresponds precisely to that value, and therefore the 
dielectric constant also needs to be evaluated at that divergent point. However our 
interest is in the radial scaling of the mobility and not on evaluating its precise value, so 
we adopt the assumption of low, but non-zero, temperature used in a similar context in 
Ref . ^1] , whereby the divergence is removed leading to a well-defined dielectric function. 
The analysis of Ref. [i4j also showed that a doubling of the wire radius caused relatively 
small changes in the dielectric constant over a wide range of carrier densities. But as 
we will establish in this paper, the effects of increasing radius on the mobility on the 
other hand is exponential in nature. Therefore the dependence of the dielectric function 
on wire size should have little qualitative impact on the radial scaling of the mobility. 
Hence we will treat the dielectric function as a constant e over the range of radii that 
we consider in this paper. 

The transverse wave functions involve Bessel functions and for a wire of radius R 
are given by 

The various channels are labelled by two indices i = l,n, with / corresponding to 
the order of the Bessel functions Ji{Kinr), and n labelling the zeroes for each order 
in a sequence of increasing magnitude. The transverse eigenmodes are determined by 
the boundary condition that the Bessel functions vanish on the surface of the wire 
Ji^KinR) = 0, with the zeroes denoted by p;„ = KinR. 

We assume a azimuthally symmetric layer of impurities of bulk density p(ro) 
distributed between radii r = a and r = b. On integrating over the impurity distribution 
we obtain the scattering amplitude sum (jSJ 

rb mL 

J2Sj,i(k\k) = L / rfro27rrop(ro)^-3— |(±W,/V|y(r,ro)|A;«„,/n)|^ 

with 

Sk'.in = ;r~ii3 [ drorop{ro) 

zne^h kiln' Ja 



Mr,0) = ^^^^^^. (11) 
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Table 1. The zeroes pin of the Bessel functions, for the lowest transverse channels 
in a uniform cylindrical wire, are listed along with the minimal values of the scaled 
wire radius n^/^R at which each channel becomes populated with charge carriers. For 
a bulk carrier density of n = 10^^ m^"^ that we use for numerical estimates, the right 
column would correspond to radii in units of 10^^ m 10 nanometers. 



P In 


n^^^R 


Po,i 


= 2.40482 





Pl.l 


= 3.83171 


0.8455 


P2.1 


= 5.13562 


0.1173 


PO.2 


= 5.52007 


0.1305 


P3.1 


= 6.38016 


0.1539 


Pi, 2 


= 7.01559 


0.1706 


P4,l 


= 7.58834 


0.1856 


P2,2 


= 8.41724 


0.2048 


P0,3 


= 8.65372 


0.2116 



The + sign corresponds to forward scattering (0 = 0) and the — to backscattering 
(0 = vr). Using this notation in Eq. (0) gives a system of hnear equations for the 
relaxation times 



{l'n'}y^{ln} 



Tin + {^l'n',ln " ^Un' ,ln) ~ (^2) 

{l'n'}j^{ln} 



At zero temperature the matrix elements, like the relaxation times, are evaluated at 
the effective Fermi wavevector kin for each channel. In the strict one dimensional limit 
when only the lowest channel is available we retrieve the well known result = 2Sqi qi 
|llj . It has been a common practice to assume that in thin wires the carriers may be 
assumed to be confined to the wire axis, thereby justifying the usage of a 5-function 
to approximate the transverse profile of the carrier density, in which case the matrix 
elements in Eq. (|T^ reduce to 

Svn'M - ^ 2t.3, / drQrop{rQ)Kl\\kin T h'n'Vo] (14) 

While this is mathematically simpler, we will presently show that this approximation is 
invalid for the cases we consider. 



4. Surface Impurity 

For some quasi-one dimensional nanowires, it is a good approximation to treat the 
impurities as distributed in a layer of varying thickness outside the wire; this models 
surface roughness or modulation doping of dopants. We therefore take the scattering 
centers to be distributed in a thin uniform layer of width w along the surface of the 
wire. First we take the bulk density of the impurities to be constant within that layer 
^('"o) = Pb- We then write the matrix elements in Eq. lfT^ in a way that makes the 
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radial dependence more transparent: 



^V'n'ln — Q ^ 



Ql'n' 



(15) 



X 



^ ^^^ Jl'{Pl'n'X)Jl{pinX) 
kJv{PVn')J'l{Pln) 



2^ de 

( 

27r 



-'^'Ko T qinWl + x^ -2xcose 



Everything that does not depend on the radius has been included in the pre-factor Q 
which contains all the dimensioned quantities. The rest of the expression contains only 
dimensionless quantities as we have rescaled the lengths by the radius R, so that the 
integration variable is x = r/R and the wavevectors are q = Rk . The scale and the 
dimension of the matrix element are then set by the constant in front 



Q 



me IqU 



e? m, 



m pb w 
n Iq 



1.04 X 10^^ s 



(16) 



1 m pbW 
z^rue n lo_ 
10~^m along the radial direction, 



We have assumed a length scale of a nanometer Iq 
for the wire radius R as well as for the width of the scattering layer w. We have also 
scaled the impurity density Pb by the carrier density n. The dielectric constant of the 
wire is denoted by . 

Another advantage of writing the matrix elements this way lies in the fact that the 
factor Q is common to all the matrix elements S^^^, in the set of linear equations in 
Eq. (|T!?jl . We can therefore divide through by that factor, and since the multiplication of 
a column of a determinant by a constant has the effect of multiplying the determinant 
by the same constant, all the relaxation times r/„ for every channel carry a common 
factor of Q^^. 

In certain cases it is more accurate to assume that the linear density of the impurities 
is constant instead of the bulk density of impurities. In that case the radial scaling is 
somewhat different; we then have to replace wpB Pl/{2ttR) so that the scattering 
matrix elements are 



c'± _ n' v (-^/^o) 

'-'I'n'ln — V ^ 



X 



q' 

^ ^^^:MP/'n'a;) Jl (pinX) 



(17) 



Jo ^e'^' ''^'^0 {\qi'n' T qmWl + x^ - 2xcos^ 



'0 2Jl'^Pl'n')JliPln) 

The coefficient Q' has the same scale factor as Q, but now the bulk impurity density 
and the width of the impurity layer are replaced by a linear density of the impurities 
pL, scaled by the bulk carrier density and a nanoscale area element n x 2t^Iq 

I m PL 



Q' 



me^lQn 



X 



rUe 2'Klln 



5. Results and Discussion 



We now proceed with numerical estimates for gallium arsenide (GaAs) for which our 
assumption of parabolicity Q is appropriate. Since we are interested mainly in the 
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radial scaling we present our results scaled by the common constant Q since it carries no 
radial dependency. But it determines the intrinsic magnitude of the physical quantities, 
therefore we first provide an estimate of its value for a GaAs nanowire. The effective 
mass for GaAs is m = 0.068me, assumed the same for all the channels, and the dielectric 
constant is about es = 12, not significantly altered by screening at the high carrier 
density that we consider ^3] . If we take the impurity density to be equal to the carrier 
density and the impurity layer to be of the order of a nanometer thick w ~ 10~^ m, 
we obtain a value of Q — 4.9 x 10^^ s^^. We may then estimate the magnitude of the 
average mobility by writing it as 



The pre-factor contains all the dimensioned quantities and, for the above-mentioned 
value of Q, it is e/{Qm) ~ 5 x 10^ cm^/V/s. Our numerical computation of 
the remaining dimensionless part, presented below, then yields a mobility in the 
range 10^ — 10'' cm^/V/s, which is consistent with the magnitudes in experimental 
measurements. 

Having established the range of magnitudes of the mobility we now turn our 
attention to the radial scaling. The numerical estimates assume a carrier density of 
10^^ cm~'^ = 10^'* m~^. When all else is fixed the radius determines the number of 
channels available for transport in Eq. The minimal radii for each channel to be 
populated with carriers having non-vanishing longitudinal energy are shown in Table I. 

In Fig. ^a) we plot the modified Bessel function of the second kind for the lowest 
channel, -K'o(x, 6'; goi) = Ko{2qQi^l + — 2x cos(6')); it contains the effects of the 
Coulomb scattering in the scattering matrix elements. In the x-9 plane, the function 
Kq[x, 9; goi) is very strongly peaked along the 9 = line, therefore we specifically chose 
to plot Ko{x, 0; goi)- In order to see the behavior over the entire range of wire radii we 
present plots for the minimum radius R = 5 nm and maximum radius R = 18.5 nm 
that we consider. Alongside we plot the un-normalized radial profile of the carrier 
density in the lowest channel Joi{Poix). We find that the overlap of the carrier density 
profile Joi{poix) with the prominent region of Ko{x,9]qoi) decreases significantly with 
increased wire radius, indicating that an approximation that replaces the radial profile 
with a delta function at the wire axis would deteriorate rapidly with larger wire radius. 
That is exactly what we see in Fig. Wljo) where we plot the intra-band scattering matrix 
element S^i ^i for the lowest channel using first Eq. (fT^ which uses the appropriate 
radial profile and secondly Eq. (|T^ which uses the ^-function approximation. The two 
curves are noticeably different even for a wire radius of 5 nm, but for larger wire radii 
they differ by several orders of magnitudes. This is exactly what one would expect for 
a surface distribution of impurities because the impurities are further removed from the 
wire axis for a larger wire. 

Figure 2 shows the scaling of the Fermi wavevector with the radius. We see that 
the dimensionless product of the wire radius with the Fermi wavevector, Rx kp has an 
almost monotonic growth with the radius. The Fermi vector itself decreases noticeably 



/i = -^9{E,) 
Qm 



Y.^q^[Qr^iqi)] 



(19) 
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Figure 1. (a) Plot of the radial profile of the carrier distribution of the lowest channel, 
Joi(poia^)i shown alongside with the radial dependence of the impurity scattering 
contained in the factors i4ro(2goi|l — x\) shown for the extreme values of the range 
of wire radii R that we consider, (b) Semi-log plot of the scattering matrix element, 
scaled by the constant Q, for the lowest channel with (thin line) and without (thick 
line) the (5-function approximation for the radial density of the carriers. 



as each new channel becomes available, and then rebounds gradually but with an overall 
decline of the peak values reached before each succeeding channel enters. So the general 
trend is that kp gets smaller with increasing radius; with sufficiently large number of 
channels we expect it to approach the bulk value which, for the carrier density we have 
assumed, would be k^P = (Stt^IO"^"^)^/^ ~ 0.31 nm~^. Our plot suggests a gradual 
approach to that limit. That limit gives a criterion for when the system makes the 
transition from quasi-lD to 3D. 

We plot the mobility in Fig. |H1 on a logarithmic scale to show that the growth 
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Figure 2. The net Fermi wavevector kp — \J 2mEF /It? ^ for the carriers in the wire 
for a fixed carrier density n = 10^'* m"^, is plotted along the left axis as a function of 
the radius. The dimensionless quantity Rkp is plotted along the right axis. 



of the mobility with radius is of an exponential nature in between points of sharp 
declines. The most striking feature is that the multi channel scattering destroys the 
simple monotonic growth of the mobility seen with a single channel [TIJ. As the radius 
increases, and each new channel becomes energetically available for scattering there is 
a sharp reduction in the average mobility. In between the addition of new scattering 
channels, the increasing wire radius causes the mobility to grow, but that growth has an 
inflection point, reflecting the competition with increased scattering. The sharp declines 
cause an overall lowering of the mobility as the radius increases significantly. This is 
consistent with an eventual approach to bulk behavior. 

We get a sense of the actual intrinsic magnitude of the mobility and its variations 
in Fig. m where a linear scale is employed. This figure also illustrates the effect of 
changing the radial scaling of the impurity density itself; we have plotted alongside the 
mobility for the case described in Eqs. (fT7|) and (fTHj) where regardless of the increased 
wire growth the linear density of the impurity layer remains constant. Those equations 
show that there is then an extra factor of R in the mobility, because the surface density 
becomes sparser with increased wire size; this causes the mobility to increase faster at 
larger radii as we see from the plot. Otherwise the two curves are quite similar implying 
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Figure 3. The mobility scaled by e/{Qm) is plotted on a logarithmic scale as a 
function of the wire radius for fixed carrier n = 10^^ m~^. The {l,n) label for each 
new channel is shown as it becomes available causing the sharp declines in the mobility. 

that the essential features are not affected significantly by the nature of the impurity 
distribution, because the strong exponential behavior of the modified Bessel functions 
Kq dominates the trend. 

We have assumed a simple model which highlights general trends in the carrier 
mobility and related parameters as a semiconductor nanowire deviates from the strict 
size quantization limit towards bulk behavior. Real semiconductor wires will not have a 
uniform cross section, and are likely to have impurities in the wire interior as well as on 
the surface, so the behavior we have seen will therefore likely be less sharp. Moreover 
it is known that screening affects intra-band scattering more significantly fH^ EBj than 
inter-band scattering. So a more accurate treatment of the dielectric function is likely to 
diminish the changes in the mobility as new channels become available. But the general 
trends we have established are sufficiently striking that those features would be relevant 
in experimentally observed behavior of the mobility as the wire size is increased. In 
particular the effect of mter-channel scattering should be manifest in the degenerate 
limit. 
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Figure 4. Plot of the scaled mobility versus the wire radius on a linear scale for 
constant bulk density of surface impurities (thick line, along the left axis) and for 
constant linear density (thin line, along the right axis). Note the scale factors are 
different for each, Q and Q' given by Eqs. and p8(l respectively. 
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